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Abstract 

We describe a parametrized Yang-Baxter equation with nonabelian param- 
eter group. That is, we show that there is an injective map g h-> R{g) from 
GL(2, C) x GL(1, C) to End(UtX'U) where V is a two-dimensional vector space 
such that if g, h G G then R\2(g)Ri3 (gh) i?23(/i) = R23{h) R13 (gh)R\2 (g) ■ Here 
Rij denotes R applied to the i, j components of V <8> V <8> V. The image of 
this map consists of matrices whose nonzero coefficients a\, 0,2, b\, 62, c\, C2 
are the Boltzmann weights for the six-vertex model, constrained to satisfy 
a\a,2 + &1&2 — c\C2 = 0. This is the exact center of the disordered regime, and 
Q\ \ is contained within the free Fermionic eight-vertex models of Fan and Wu. As 

an application, we give a new proof based on the Yang-Baxter equation of a 
result of Hamel and King representing a Schur polynomial times a deforma- 
tion of the Weyl denominator as the partition function of a six-vertex model. 
$_i ■ Furthermore, the parameter group can be expanded (within the eight-vertex 



model) to a group having GL(2) x GL(1) as a subgroup of index two. In this 
expanded context we find a second representation of Schur polynomials times 
a different deformation of the Weyl denominator as a partition function. These 
structures give a Yang-Baxter system in the sense of Hlavaty. 

Baxter's method of solving lattice models in statistical mechanics is based on the 
star-triangle relation, which is the identity 

-^12*513^23 = T23S13-R12, (1) 
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where R, S, T are endomorphisms of V <E> V for some vector space V. Here Rij is the 
endomorphism of V <E> <E> V in which i? is applied to the z-th and j-th copies of V 
and the identity map to the k-th component, where i,j,k are 1,2,3 in some order. If 
the endomorphisms R,S,T are all equal, this is the Yang-Baxter equation (cf. [PT] . 

E3). 

More generally, one may ask for solutions to a parametrized Yang-Baxter equa- 
tion, where the endomorphism R now depends on a parameter g (ranging over a 
group G) and (CQ) takes the form 

RMRu(9 ■ h)R 2 s(h) = R 23 (h)R 13 (g ■ h) R 12 (g) (2) 

for arbitrary choice of parameters g,h G G. There are many such examples in the 
literature in which the group G is an abelian group such as R or R x . In this paper 
we present an example of (j2J) having a non-abelian parameter group. The example 
arises from two-dimensional lattice models - the six- and eight-vertex models. 

We now briefly review the connection between lattice models and instances of 
(Op) and (j2j). In statistical mechanics, one attempts to understand global behavior 
of a system from local interactions. To this end, one defines the partition function 
of a model to be the sum of certain locally determined Boltzmann weights over all 
admissible states of the system. Baxter (see [1] and [2], Chapter 9) recognized that 
instances of the star-triangle relation allowed one to explicitly determine the partition 
function of a lattice model. 

The six-vertex, or 'ice-type,' model is one such example that is much studied 
in the literature, and we revisit it in detail in the next section. For the moment, 
we offer a few general remarks needed to describe our results. In our presentation 
of the six-vertex model, each state is represented by a labeling of the edges of a 
finite rectangular lattice by ± signs, called spins. If the Boltzmann weights are 
invariant under sign reversal the system is called field-free, corresponding to the 
physical assumption of the absence of an external field. For field-free weights, the 
six- vertex model was solved by Lieb [25] and Sutherland |36j . meaning that the 
partition function can be exactly computed. The papers of Lieb, Sutherland and 
Baxter assume periodic boundary conditions, but non-periodic boundary conditions 
were treated by Korepin [20] and Izergin [16] . Much of the literature assumes that the 
model is field-free. In this case, Baxter shows there is one such parametrized Yang- 
Baxter equation with parameter group C x for each value of a certain real invariant 
A, defined below in (Q in terms of the Boltzmann weights. 

One may ask whether the parameter subgroup C x may be enlarged by including 
endomorphisms whose associated Boltzmann weights lie outside the field-free case. 
If A 7^ the group may not be so enlarged. However we will show in Theorem [3] that 
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if A = 0, then the group C x may be enlarged to GL(2, <D) x GL(1, C) by expanding 
the set of endomorphisms to include non-field-free ones. In this expanded A = 
regime, R{g) is not field-free for general g. It is contained within the set of exactly 
solvable eight- vertex models called the free Fermionic model by Fan and Wu [6], [7]. 
Our calculations suggest that it is not possible to enlarge the group G to the entire 
free Fermionic domain in the eight vertex model but we are able to enlarge G to a 
group containing GL(2, C) x GL(1, C) as a subgroup index two (Theorem E]). 

In Section |2] we give a heuristic argument to show that if there is a set of 
endomorphisms such that for any S and T in that set there exists R such that 
R12S13T23 = T23S13R12 then an associativity property is satisfied, so that fl2J) is satis- 
fied. Of course our rigorous results do not depend on this plausible reasoning, but it 
seems useful to know that the associativity that we observe is not entirely accidental. 

As an application of these results, we study the partition function for ice-type 
models having boundary conditions determined by an integer partition A and Boltz- 
mann weights chosen so that both A = and so that the degenerate case A = 
matches the standard deformation of Weyl's denominator formula for GL n (<D). This 
leads to an alternate proof of a deformation of the Weyl character formula for GL n 
found by Hamel and King [12], [TT]. That result was a substantial generalization 
of an earlier generating function identity found by Tokuyama [37], expressed in the 
language of Gelfand-Tsetlin patterns. 

More precisely, we will exhibit two particular choices of Boltzmann weights and 
boundary conditions in the six-vertex model giving systems and &^ for every 
partition A of length ^ n. We will prove that the partition functions are 

z {&\) = J\(Uzj + Zi)sx{zi, ■■■ ,z n ), z {&\) = Y\{tjZj + Zi)sx{zi, ■■■ ,z n ), (3) 

i<j i<j 

where ti are deformation parameters and s\ is the Schur polynomial (Macdonald [26]). 
The method of proof is inspired by ideas of Baxter in [1] and [2], though the Boltz- 
mann weights we use are not field-free. The A model is essentially that given by 
Hamel and King. The notation here is somewhat unfortunate as A denotes a recipe 
for choosing weights and A denotes an invariant defined in terms of weights, but has 
been chosen to match earlier uses of this notation in the literature. 
To justify these evaluations of the partition function define 

S\ [Zi, ■■■ ,z n ;ti, ••• ,t n ) — = — 77— — — r- ( 4 ) 

LLi<j\ Z i Z j 1" Z i> 

Then one seeks to show that s\ is symmetric in the sense that it is unchanged if the 
same permutation is applied to both and t^. Once this is known, it is possible to 
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show that it is a polynomial in the Z{ and U, then that it is independent of the U; 
finally, taking ti — — 1 one may invoke the Weyl character formula and conclude that 
it is equal to the Schur polynomial. 

In order to prove the symmetry property of we will use an instance of fl2]) with 
A = 0. We thus obtain a new proof of Tokuyama's formula and of Corollary 5.1 
in Hamel and King [12] , which is our Theorem [TTJ A second instance of the star- 
triangle relation solves the same problem for the analogously defined s^, and a third 
instance shows directly without using the above evaluations, that s\ = . 

There are, as we have mentioned, Boltzmann weights of two different types V and 
A. (We refer to these as different types of "ice.") Moreover if X, Y G {T, A} we will 
give an R-matrix Rxy which has the effect of interchanging a strand of X ice with a 
strand of Y ice; thus in (TTJ , S is of type X and T is of type Y. We will prove that the 
R-matrices Rrr and Raa both satisfy the Yang-Baxter equation, and we will prove 
similar relations that involve all four types of ice Rxy i n various combinations. 

Of the six types of ice that we will consider: T, A, Rrr, Rta, Raf and Raa, 
only T and Rrr come from the space of endomorphisms parametrized by GL(2, C) x 
GL(1, C). The others may be accommodated by enlarging the parameter group to a 
disconnected group having GL(2, <C) x GL(1, <D) as a subgroup of index two. 

In another direction, Hlavaty [13] has defined the notion of a Yang-Baxter system. 
As in our setup, this involves six types of endomorphisms. His definition has two 
independent motivations. On the one hand, there is the work of Freidel and Mail- 
let [10] on integrable systems, and on the other hand, there is work of Vladimirov [39] 
which attempts to clarify the relation of the construction of Faddeev, Reshetikhin 
and Takhtajan [5] to Drinfeld's quantum double. In Section [0J we show that our 
construction is an example of a Yang-Baxter system. In the case where the ti are 
equal, these Yang-Baxter systems are related to those previously found by Nichita 
and Parashar [31], [50] . 

Our boundary conditions depend on the choice of a partition A. Once this choice 
is made, the states of the model are in bijection with strict Gelfand-Tsetlin patterns 
having a fixed top row. These are triangular arrays of integers with strictly decreasing 
rows that interleave (Section H]). Since in its original form Tokuyama's formula 
expresses what we have denoted Z{G T X ) as a sum over strict Gelfand-Tsetlin patterns 
it may be expressed as the evaluation of a partition function. 

This connection between states of the ice model and strict Gelfand-Tsetlin pat- 
terns has one historical origin in the literature for alternating sign matrices. (An 
independent historical origin is in the Bethe Ansatz. See Baxter [2] Chapter 8 and 
Kirillov and Reshetikhin [TO].) The bijection between the set of alternating sign 
matrices and strict Gelfand-Tsetlin patterns having smallest possible top row is in 
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Mills, Robbins and Rumsey [29], while the connection with what are recognizably 
states of the six- vertex model is in Robbins and Rumsey [33]. This connection was 
used by Kuperberg [21] who gave a second proof (after the purely combinatorial one 
by Zeilberger [30]) of the alternating sign matrix conjecture of Mills, Robbins and 
Rumsey [29]. Kuperberg' s paper follows Korepin [20J and Izergin [16] and makes use 
of the Yang-Baxter equation. It was observed by Okada [32] and Stroganov [35J that 
the number of n x n alternating sign matrices, that is, the value of Kuperberg's ice 
(with particular Boltzmann weights involving cube roots of unity) is a special value 
of the particular Schur function in 2n variables with A = (n, n, n — 1, n — 1, • • • ,1,1) 
divided by a power of 3. Moreover Stroganov gave a proof using the Yang-Baxter 
equation. This occurrence of Schur polynomials in the six-vertex model is different 
from the one we discuss, since Baxter's parameter A is nonzero for these investiga- 
tions. 

There are other works relating symmetric function theory to vertex models or 
spin chains. Lascoux [23J, [22] gave six- vertex model representations of Schubert 
and Grothendieck polynomials of Lascoux and Schiitzenberger [21] and related these 
to the Yang-Baxter equation. Fomin and Kirillov [8], [9] also gave theories of the 
Schubert and Grothendieck polynomials based on the Yang-Baxter equation. Tsile- 
vich [38] gives an interpretation of Schur polynomials and Hall-Littlewood polyno- 
mials in terms of a quantum mechanical system. Jimbo and Miwa [18] give an 
interpretation of Schur polynomials in terms of two-dimensional Fermionic systems. 
(See also Zinn- Justin [4T].) 

McNamara [28J has clarified that the Lascoux papers are potentially related to 
ours at least in that the Boltzmann weights [23] belong to the expanded A = 
regime. Moreover, he is able to show based on Lascoux' work how to construct 
models of the factorial Schur functions of Biedenharn and Louck. 

We are grateful to Gautam Chinta and Tony Licata for stimulating discussions. 
This work was supported by NSF grants DMS-0652609, DMS-0652817, DMS-0652529 
and DMS-0702438. SAGE [34J was very useful in the preparation of this paper. 

1 The Six- Vertex Model 

We review the six-vertex model from statistical mechanics. Let us consider a lattice 
(or sometimes more general graph) in which the edges are labeled with "spins" ±. 
Depending on the spins on its adjacent edges, each vertex will be assigned a Boltz- 
mann weight. 

The Boltzmann weight will be zero unless the number of adjacent edges la- 
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beled '— ' is even. Let us denote the possibly nonzero Boltzmann weights as follows: 
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We will consider the vertices in two possible orientations, as shown above, and ar- 
range these Boltzmann weights into a matrix as follows: 



R 



( oi di\ 
h ci 
c 2 b 2 
\d 2 a 2 J 



( ai(R) di(R) \ 

h(R) a{R) 

c 2 {R) h{R) 
V d 2 (R) a 2 (R) J 



(5) 



If the edge spins are labeled v, (3, 7, 9 £ {+, — } as follows: 




then we will denote by R \ the corresponding Boltzmann weight. Thus R\X = a i(-^); 
etc. Because we will sometimes use several different systems of Boltzmann weights 
within a single lattice, we label each vertex with the corresponding matrix from 
which the weights are taken. 

Alternately, R may be thought of as an endomorphism of V <g> V, where V is a 
two-dimensional vector space with basis v + and t>_. Write 



v+, v. 



(6) 

v_ gives ()6]) as 



Then the ordering of basis vectors: v + (g) u + , u + £g> u_, w_ 
the matrix (JSJ). 

If is an endomorphism of V <S> V we will denote by 0i 2 , 0i3 and 023 endomor- 
phisms oiV ®V ®V defined as follows. If = 0' <g) 0" where 0', 0" £ End(K) then 
1, 0i3 = 0' (g> 1 <g> 0" and 023 = 1 (g) 0' <E> 0". We extend this definition 
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to all (p by linearity. Now if <j),ip,x are three endomorphisms of V £g> V we define the 
Yang-Baxter commutator 

[0, ^, Xl = 012^13X23 - X23^13012- 

Lemma 1 The vanishing of [i?, 5, T] is equivalent to the star-triangle identity 




for every fixed combination of spins a, r, a, (3, p, 6. 

The term star-triangle identity was used by Baxter. The meaning of equation 
(0) is as follows. For fixed a, r, a, /3, p, 9, //, u, 7, the value or Boltzmann weight of 
the left-hand side is just the product of the Boltzmann weights at the three vertices, 
that is, R^SJpTftj, and similarly the right-hand side. Hence the meaning of (J7J) is 
that for fixed a, r, a, (3, p, 6, 

£ rv £S%t% = £ TiSs£n>&. (8) 

Proof Let us apply {R, S, T] to the vector v a ® v T ® Vp. On the one hand by © 

^12^13723(^0-®^®^) = ^12^13 ^ (Vg ®Vj,® Vs) 

ip,8 

= R12 £ s ip T is( v <i>® v i>® v <x) 

ip,8,(f>,a 
il>,8,<f>,a,6,p 

and similarly 

V0 ®v p ® v a ). 
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We see that the vanishing of {R, S, T] is equivalent to (JHI) • □ 

In this section we will be concerned with the six-vertex model in which the weights 
are chosen so that d\ = d 2 = in the table above. In [2], Chapter 9, Baxter considered 
conditions for which, given S and T, there exists a matrix R such that [i?, S, T] = 0. 
We will slightly generalize his analysis. He considered mainly the field-free case 
where en (72) = a 2 (R) = a(R), h(R) = b 2 (R) = b(R) and d(R) = c 2 (R) = c(R). The 
condition C\(R) = c 2 (R) = c(R) is easily removed, but with no gain in generality. 
The other two conditions a\{R) = a 2 (R) = a(R), bi(R) = b 2 (R) = b(R) are more 
serious restrictions. 

In the field-free case, let 

A(i^) = Q( ^ )2 2 t ( ^(^) C(jR)2 - ai(R) = a 2 (R) = a(R), etc. (9) 

Then Baxter showed that given any S and T with A(S') = A(T), there exists an R 
such that {R,S,Tj = 0. 

Generalizing this result to the non-field-free case, we find that there are not one 
but two parameters 

a 1 (R)a 2 (R) + h(R)b 2 (R) - Cl (R)c 2 (R) 



MR) 
MR) 



2a 1 (R)b 1 (R) 
a 1 (R)a 2 (R) + b^R^R) - Cl (R)c 2 (R) 



2a 2 (R)b 2 (R) 
to be considered. 

Theorem 1 Assume that ai(S), a 2 (S), bi(S), b 2 (S), ci(S), c 2 (S), ai(T), a 2 (T), 
b\(T), b 2 (T), Ci(T) and c 2 (T) are nonzero. Then a necessary and sufficient condition 
for there to exist parameters a\(R), a 2 (R), b\(R), b 2 (R), c\(R), c 2 (R) such that 
{R,S,Tj = withc 1 {R),c 2 {R) nonzero is that A^S) = Ai(T) and A 2 (S) = A 2 (T). 

Proof Suppose that A 1 (S') = Ai(T) and A 2 (S) = A 2 (T). Then we may take 
ai (R) 

ai{i ) 



a 2 (R) 



b 2 (S)ai(T)bi(T) ■ 


- ai(S)bi(T)b 2 (T) + ai (S) Cl (T)c 2 (T) 


ai (505i(5)o2(T) ■ 


oi(T) 

- a 1 {S)a 2 {S)b 1 {T) + c 1 {S)c 2 {S)b 1 {T) 




h(S) 


h(S)a 2 (T)b 2 (T) 


- a 2 (S)h{T)b 2 (T) + a 2 (S) Cl {T)c 2 (T) 


a 2 (S)b 2 (S) ai (T) 


a 2 (T) 

- ai {S)a 2 (S)b 2 (T) + cx(S)c 2 (S)b 2 (T) 



b 2 (S) 
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h(R) = h(S)a 2 (T) - a 2 {S)bi(T), b 2 (R) = b 2 (S)a 1 (T) - ai(S)b 2 (T), (12) 

Cl (R) = Cl (S)c 2 (T), c 2 (R) = c 2 (S) Cl (T). (13) 

Using Ai(S') = Ai(T) and A 2 (S) = A 2 (T) it is easy to that the two expressions 
for ai(R) agree, and similarly for a 2 (R). One may check that [i?, S, T] = 0. On 
the other hand, it may be checked that the relations required by {R, S,TJ =0 are 
contradictory unless Ai(S') = Ai(T) and A 2 (<5') = A 2 (T). □ 

In the field- free case, these two relations reduce to a single one, A(S') = A(T), 
and it is remarkable that A(i?) has the same value: A(i?) = A(S') = A(T). 

This equality has important implications for the study of row-transfer matrices, 
one of Baxter's original motivations for introducing the star-triangle relation. Given 
Boltzmann weights ai(R), a 2 (R), ■ ■ ■ , we associate a 2" x 2" matrix V(R). The 
entries in this matrix are indexed by pairs a = (a%, • • - , a n ), /?=(/?!,■••, (3 n ), where 
cti,(3i G {±}. If G {±} we may consider the Boltzmann weight of the 

configuration: 



a.\ ct 2 a 



£l « S 2 ■ 



Here e n+ i = £i, so the boundary conditions are periodic. The coefficient V(R) a ,p is 
then the "partition function" for this one-row configuration, that is, the sum over 
possible states (assignments of the 

It follows from Baxter's argument that if R can be found such that [R, S,TJ = 
then V(S) and V(T) commute, and can be simultaneously diagonalized. We will 
not review Baxter's argument here, but variants of it with non-periodic boundary 
conditions will appear later in this paper. 

In the field-free case when l_R, S, T] = 0, V(R) belongs to the same commuting 
family as V(S) and V(T). This gives a great simplification of the analysis in Chap- 
ter 9 of Baxter [2] over the analysis in Chapter 8 using different methods based on 
the Bethe Ansatz. 

In the non-field-free case, however, the situation is different. If Ai(S') = Ai(T) 
and A 2 (S) = A 2 (T) then by Theorem [1] there exists R such that {R, S, T] = 0, 
and so one may use Baxter's method to prove the commutativity of V(S) and V(T). 
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However Ai(R) and A 2 (R) are not necessarily the same as Ai(S') = Ax(T) and 
A 2 (S) = A 2 (T), respectively, and so V(R) may not commute with V(S) and V{T). 

In addition to the field-free case, however, there is another case where V(R) 
necessarily does commute with V(S) and V(T), and it is that case which we turn to 
next. This is the case where aia 2 + b\b 2 — Cic 2 — 0. The next theorem will show that 
if the weights of S and T satisfy this condition, then R exists such that {R, S,TJ = 0, 
and moreover the weights of R also satisfy the same condition. Thus not only V(S) 
and V(T) but also V(R) lie in the same space of commuting transfer matrices. 

In this case, with a\ = ai(R), etc., we define 



n(R) 



TV 



bi ci 
c 2 b 2 



\ 



a 2 J 



l>2 

a 2 



\ 



c 2 J 



(14) 



Theorem 2 Suppose that 



ai(S)o2(S) + hiSMS) - Cl (S)c 2 {S) = ai {T)a 2 (T) + h{T)b 2 {T) - Cl (T)c 2 (T) = 0. 

(15) 

Then the R G End(V <g> V) defined by n(R) = ^(^^(r)- 1 satisfies {R,S,T} = 0. 
Moreover, 

ai(R) a 2 (R) + h(R) b 2 (R) - a(R) c 2 (R) = 0. (16) 

Proof The matrix R will not be the matrix in Theorem [H but will rather be a 
constant multiple of it. We have 



vrm- 1 = - 

V ; D 



( c 2 (T) 



\ 



a 2 {T) -b 2 {T) 
h(T) ai (T) 



\ 



Ci(T) J 



where D = ai(T)a 2 (T) + bi(T)b 2 (T) = ci(T)c 2 (T). With notation as in Theorem [TJ 
using ( TH)]) equations (1TU1) and ffTTl) may be written 

ai (R) = a 1 (S)a 2 (T) + b 2 (S)b 1 (T), 
a 2 (R) = a^a^ + b^SMT). 

Combined with (Tl2l) and ffl3l) these imply that tt(R) = vr(S') Dir(T)~ l . However we 
are free to multiply R by a constant without changing the validity of \R, S,TJ = 0, 
so we divide it by D. □ 
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We started with S and T and produced R such that {R, S,TJ = because this 
is the construction motivated by Baxter's method of proving that transfer matrices 
commute. However it is perhaps more elegant to start with R and T and produce S 
as a function of these. Thus let 7Z be the set of endomorphisms R of V <8> V of the 
form (jSJ) where aia 2 + bib 2 = Cic 2 . Let 72* be the subset consisting of such R such 
that c\Ci ^ 0. 

Theorem 3 There exists a composition law on 1Z* such that if R,T 6 1Z* , and if 
S = RoT is the composition then [i?, T] =0. 27ms composition law is determined 
by the condition that ir(S) = tt(R)tt(T) where n : 1Z* — > GL(4, C) is the map 
Then 1Z* is a group, isomorphic to GL(2, C) x GL(1, C). 

Proof This is a formal consequence of Theorem [2j □ 

It is interesting that, in the non-field-free case, the group law occurs when 
Ai = A 2 = 0. In the application to statistical physics for field-free weights, phase 
transitions occur when A = ±1. If |A| > 1 the system is "frozen" in the sense 
that there are correlations between distant vertices. By contrast — 1 < A < 1 is the 
disordered range where no such correlations occur, so our group law occurs in the 
analog of the middle of the disordered range. 

2 Composition of R-matrices 

Theorem [31 defining a group structure on a set of R-matrices, may be regarded as a 
non-abelian parametrized Yang-Baxter equation. In our example, the composition 
law on R-matrices that makes S the product of R and T when [i?, S, T] =0 is 
associative because of its definition in terms of matrix multiplication. In this section, 
we give a heuristic argument suggesting that any time we have such a composition 
law defined by the vanishing of a Yang-Baxter commutator, associativity should 
follow. This section is not needed for the sequel. 

Let us assume that we are given a vector space V over a field F and a subset 7Z 
of End(V ®V) which is homogeneous in the sense that if ^ R G 1Z then 1Z contains 
the entire ray FR. Let P(JZ) be the set of such rays. 

Let us assume that if R, T are nonzero elements of 7Z then there is another S € 7Z 
that is unique up to scalar multiple such that {R, S, T] = 0. As we remarked before 
Theorem [2] there might be such an S that would be useable for applications but 
that it might not lie in the same space 7Z, and indeed this is the usual situation 
for the six vertex model with weights that are not field-free and also not in the free 
Fermionic case of Theorem [2J But with this assumption, (R, T) i— > S is a well-defined 
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composition law on P(1Z). Let us denote this composition S = R o T. We will give 
a plausible argument that this composition law should be associative. 

We begin with three nonzero elements R, S, T of 1Z. We will compare endomor- 
phisms of V <8> V <S> V <S> V. In addition to identities such as R\ 2 {R ° S) 13823 = 
S 23 (R o S)i 3 Ri2 we will use identities such as R\ 3 T 24 = T 24 R 13 which are true for 
arbitrary endomorphisms of V <S> V. Let 

X 1U = (R o S) 13 (R o (S o T)) 14 T 34 , X( 34 = T 34 (R o (S o T)) 14 (R o S) 13 . 

First, we have S 2 3(S o T) 24 Xi 34 ii!i 2 equal to 

S 23 (S T) 24 (R o S^i? o(5o T))i4T 34 i?i2 = 

S 23 (R o S) 13 (S o T) 24 (R o{S o T)) u R 12 T u = 

S 23 (R o 5)i 3 i2i 2 (i2 °(So T)) 14 (S o T) 24 T 34 = 

i?i 2 (-R o S) 13 S 23 (R o(So T))u(S o T) 24 T 34 = 

#12 (-R o S^i? o(S o T)) 14 S 23 (S o T) 24 T 34 = 

R 12 (Ro S) 13 (Ro (S oT)) 14 T 34 (S oT) 24 S 23 = Ri 2 X lu (S o T) 24 S 23 . 

Using another string of manipulations, we have S 23 (S o T) 24 X' 1M R 12 equal to 

S 23 (S o T) 24 T 34 (R o (S o T)) 14 (R o 5) 13J R 12 = 
T 34 (S o T) 24 S 23 (R o (5 o T)) 14 (i? o 5) 13 i2 12 = 
T 34 (S o T) 24 (i? o (5 o T)) U S 23 (R o 5) 13J R 12 = 
T 34 (S o T) 24 (i? o (5 o T)) 14 R 12 (R o 5) 13 5 23 = 
T 34 R l2 (R o (5 o T)) 14 (S o T) 24 (R o 5) 13 5 23 = 

Ri2T 34 (Ro(SoT)) 14 (RoS) 13 (SoT) 24 S 23 = R 12 X' 1U (S o T) 24 S 23 . 

Now consider an endomorphism X of V ® V ® V that is constrained to satisfy 

S 23 (S o T) 24 Xi 34 R 12 = Ri 2 X 134 (S o T) 24 S 23 . 

This is a linear equation in the matrix coefficients of X in which the number of 
conditions exceeds the number of variables. It is reasonable to assume that if this 
has a nonzero solution that solution is determined up to constant multiple. Therefore 
(up to a constant) we have X\ 34 = X' 1U , that is, 

(RoS) 13 (Ro(SoT)) 

14^34 = T 34 (R o(So T)) 14 (R o 5)13. 

Taking the determinant shows that the constant must be a root of unity. If F — R 
or C and P(1Z) is connected, then by continuity this constant must be 1. This 
means that (R o (S o T)) satisfies the definition of (Ro S) o T, so at least plausibly, 
a composition law defined this way should be associative. 
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3 Gamma ice 

Let zi, • • • ,z n and ti, ■ • • ,t n be complex numbers with all % 7^ 0. We will refer to 
the Zi as spectral parameters and the ti as deformation parameters since these are 
the roles these variables will play when we turn to Tokuyama's theorem. Denote 



rd) 



( i 



V 



U Ziiti + 1) 
1 Zi 



\ 



Zi J 



TrrU 



/ Zi{ti + 1) \ 
1 Zi 

V 1 / 



Let 7Trr{i,j) = const X7Tr(i)vrr(i) 1 , where it is convenient to take the constant to 
be Zj(tj + 1). It follows from Theorem [2] that 

[#rr(M'),r(z),r(j)] = 0, (17) 

where Rrr(i,j) is related to 7rrr(^, j) by the relation (fl4l) . Concretely, 

\ 

(18) 

Zi ~\~ tiZj J 

The six types of vertices corresponding to the non-zero entries of T(i) and R-prihj) 
are given in Table [IJ together with their Boltzmann weights. 

Theorem 4 The star-triangle identity 



R 



IT 



' Zj ~\~ tj Z% 



V 



t^Zj Zi\p% ~t~ l y 

Zj ^ t j I X ^ z^ z j 



y 




E 



S,cf>,ip 




is valid with Boltzmann weights as in TableUl 

Proof This follows from Theorem [2] since 'n'rrihj) = const X7Tr(i)7rr(j) _1 . □ 
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• e^ ^ 

i .7 
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9 ~Q\ XD • 
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9 ~®\ /O • 
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■ UA AST 9 


• 0\ /O • 


Boltzmann 
weight 


2^ — h 2<j 




tiZj tjZi 


2^ 


(*i + 1)* 


+ l)z j 



Table 1: Boltzmann weights for Gamma ice and Gamma-Gamma ice. 

We will use Gamma ice to represent Schur polynomials, which are essentially 
the characters of finite-dimensional irreducible representations of GL n (C). If p = 
(fjLi, • • ■ , p n ) G J/ 1 then we may regard p as an element of the GL n (C) weight lattice 
and call it a weight. li p 1 ^ ■ ■ • ^ p n we say it is dominant, and if p 1 > ■ ■ • > p n 
we say it is strictly dominant. If p is dominant and p n ^ 0, it is a partition. 

Note: The word "partition" occurs in two different senses in this paper. The 
partition function in statistical physics is different from partitions in the combina- 
torial sense. So for us a reference to a "partition" without "function" refers to an 
integer partition. Also potentially ambiguous is the term "weight," referring to an 
element of the GL n weight lattice, which we identify with 7L n . Therefore if we mean 
Boltzmann weight, we will not omit "Boltzmann." 

Let A = (Ai, • • • , A n ) be a fixed partition. We will denote p — (n — l,n — 2, • • • , 0). 
We will consider a rectangular grid with n rows and \\+n columns. We will number 
the columns of the lattice in descending order from Ai + n — 1 to 0. 

A state of the model will consist of an assignment of "spins" ± to every edge. 
We will also assign labels to the vertices themselves, which will be integers between 
1 and n. For Gamma ice the vertices in the z-th row will have the label i. The spins 
of the boundary edges are prescribed as follows. 

Boundary Conditions determined by A. On the left and bottom boundary edges, 
we put +; on the right edges we put —. On the top, we put — at every column labeled 
Aj + n — % (1 ^ % ^ n), that is, for the columns labeled with values in A + p. Top 
edges not labeled by A« + n — i for any i are given spin +. 
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For example, suppose that n = 3 and A = (3, 1, 0), so that A + p = (5,2,0). Then 
the spins on the boundary are as in the following figure. 



5 4 3 2 1 


e-fo •<> - wo- - we 



K>K>K>K>K>fe 



3 o 



3 o 



3 o 



© © © © © © 



3 o 



3 o 



3 © 



(19) 



The column labels are written at the top, and the vertex labels are written next to 
each vertex. The edge spins are marked inside circles. We have left the edge spins on 
the interior of the domain blank, since the boundary conditions only prescribe the 
spins we have written. The interior spins are not entirely arbitrary, since we require 
that at every vertex "•" the configuration of spins adjacent to the vertex be one of 
the six listed in Table [1] under "Gamma ice." 

Let (5^ be the Gamma ensemble determined by A, by which we mean the set of all 
such configurations, with the prescribed boundary conditions. If x G &\, we assign 
a value w(x) called the Boltzmann weight. Indeed, Table Q] assigns a Boltzmann 
weight to every vertex, and w(x) is just the product over all the vertices of these 
Boltzmann weights. The partition function Z(&) of an ensemble (3 is $^ z6S w i x )- 
As an example, suppose that n = 2 and I = (0, 0) so A + p = (1, 0). In this case (3^ 
has cardinality two, and Z(&\) = t\Zi + Z\. The states are: 



state 


©fofb 
© © 


Q Q 
©^ogo 

©1©jlO 

© © 


Boltzmann weight 


t\Z 2 





The partition function for general A of arbitrary rank will be evaluated later in this 
paper using the star-triangle relation. 
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4 Gelfand-Tsetlin patterns 



Let us momentarily consider a Gamma ice with just one layer of vertices, so there 
are three rows of spins. Let a%, ■ ■ ■ , a m be the column numbers (from left to right) 
of — 's in the top row of spins, and let • • • , P m ' be the column numbers of — 's in 
the bottom row of edges. For example, in the ice 

5 4 3 2 1 

Q $ $ Q $ Q 
(D « « « (D « « « 

6 6 6 6 6 6 

we have m = 3, m! = 2, (ax,a 2 ,a 3 ) = (5,2,0) and (fix, fi 2 ) = (3,0). Since the 
columns are labeled in decreasing order, we have ol\ > a 2 > • ■ ■ an d /3i > (3 2 > ■ ■ ■ . 

Lemma 2 Suppose that the spin at the left edge is +. Then we have m = m! or 
m' + 1 and a\ ^ /3i ^ a 2 ^ . . . . If m = m' then the spin at the right edge is +, 
while if m — m' + 1 it is —. 

We express the condition that aj ^ /?i ^ a 2 ^ • • • by saying that the sequences 
ax, a 2 , ■ ■ ■ and fix, (3 2 , ■ ■ ■ interleave. This Lemma is essentially the line-conservation 
principle in Baxter [2], Section 8.3. 

Proof The spins in the middle row are determined by those in the top and bottom 
rows and the left-most spin in the middle row, which is +, since the edges at each 
vertex have an even number of + spins. If the rows do not interleave then one of the 
illegal configurations 



e 
e 



e 
e-«-e 

e 



will occur. Thus «i ^ (3x since if not, the vertex in the fix column would be sur- 
rounded by spins in the first illegal configuration. Now fix ^ ct 2 since otherwise the 
vertex in the a 2 column would be surrounded by spins in the second above illegal 
configuration, and so forth. The last statement is a consequence of the observation 
that the total number of spins must be even. □ 

We recall that a Gelfand-Tsetlin pattern is a triangular array of dominant weights, 
in which each row has length one less than the one above it, and the rows interleave. 
The pattern is called strict if the rows are strictly dominant. 
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It follows from Lemma [2] that taking the locations of — in the rows of vertical 
lattice edges gives a sequence of strictly dominant weights forming a strict Gelfand- 
Tsetlin pattern. For example, given the state 



ffi ffi 

0)^ToT(DioToTo 

© 9 © 9 

o^o^o^o^o^oR) 

O, 0, 0, 0o 0, 0, 

e e e e e e 

the corresponding pattern is 




(20) 



It is not hard to see that this gives a bijection between strict Gelfand-Tsetlin patterns 
and states with boundary conditions determined by A. Let us say that the weight of 
a state is (fj,i, ■ ■ ■ , /j, n ) if the Boltzmann weight is the monomial z M = n zf* times a 
polynomial in If % is a Gelfand-Tsetlin pattern, let c4(T) be the sum of the fc-th 
row. We let d n+1 (%) = 0. 

Lemma 3 If % is the Gelfand-Tsetlin pattern corresponding to a state of weight \x, 
then fj,k = d k (Z) - d k+1 (%). 

Proof From Table HJ \i k is the number of vertices in the fc-th row that have an edge 
configuration of one of the three forms: 





e^e 




6 


6 


6 



Let cti's (respectively /3j's) be the column numbers for which the top edge spin (respec- 
tively, the bottom edge spin) of vertices in the k-th row is — (with columns numbered 
in descending order, as always). By Lemma [2] we have «i ^ j3\ ^ «2 ^ • • • ^ ctn+x-k- 
It is easy to see that the vertex in the j-column has one of the above configurations 
if and only if its column number j satisfies aij > j ^ for some i. Therefore the 
number of such j is X) «» - X) A = ^Jfe(^) — ^fc+iO^)- ^ 
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5 Evaluation of Gamma Ice 



In this section we will prove the following result. 
Theorem 5 Let A = (Ai, • • • , A n ) be a partition. Then 



To begin with, define 



i<j 



(21) 



We will eventually show that is the Schur polynomial s\. But a priori it is not 
obvious from this definition that is symmetric, nor that it is a polynomial, nor 
that it is independent of t. 

Lemma 4 The expression (tk+iZk + Z]- + i)Z(&\) is invariant under the interchange 
of the spectral and deformation parameters: (zk,tk) < — > (^+1,^+1)- 

Proof We modify the ice by adding a Gamma-Gamma R-vertex (that is, one of the 
vertices from the bottom row in Table [1]) to the left of the k and k + 1 rows. Thus 
(|T9j) becomes (with k = 2 for illustrative purposes) 

5 4 3 2 1 

e e e e e e 
eSo2cA>2o2o2e 




cts ft\ rf\ ff\ rf\ ff\ 

VP \U vD \U vD vD 



which is a new boundary value problem. The only legal values for a and b are +, so 
every state of this problem determines a unique state of the original problem, and 
the partition function for this state is the original partition function multiplied by 
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the Boltzmann weight of the R- vertex, which is tk+iZk + z k +i- Now we apply the 
star-triangle identity, and obtain equality with the the following configuration. 

5 4 3 2 1 

e 00000 
e3-<^cA>2o2o2e 

rAr^ rArArArArA^ 




A\ ff\ ff\ ff\ 

VP \D \D \D \U 



Thus if 6' denotes this ensemble the partition function Z(&) = (t k+ iz k +z k+ \)Z(G\). 

Repeatedly applying the star-triangle identity, we eventually obtain the configu- 
ration in which the R- vertex is moved entirely to the right. 

5 4 3 2 1 

000000 

clclclclclq 




sfcs As A\ A\ A\ A\ 
\U \U \U \U \& \D 



Now there is only one legal configuration for the R- vertex, so c = d — — . The 
Boltzmann weight at the R-vertex is therefore t k z k+ \ + z k . Note that (z k ,tk) and 
(z k+1 ,t k+1 ) have been interchanged. This proves that (t k+1 z k + z k+ i)Z{& T x ) is un- 
changed by switching (z k ,t k ) and (z k+1 ,t k+1 ). □ 



Proposition 1 s\ is a symmetric polynomial in z\, 
the ti. 



z n , and is independent of 



Proof Consider 



Wit^ + z^Zi&l). 



(22) 



i<j 



We will show that this is invariant under the interchange k •H- k + 1. This means 
that we interchange both z k with z k+ i and t k with t k+ i. Indeed, we may write ( 122|) 
as (t k+ \z k + z k+ i)Z(&\) times the product of all factors tjZi + Zj with i < j except 
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(i, j) = (k, k + 1). These factors are permuted by k <H- k + 1, so the statement follows 
from Lemma HJ Thus ( 1221) is invariant under permutations of the indices, where it 
is understood that the same permutation is applied to the tj as to the z%. Now ( 1221) 
equals Ylijtj(^j z i + s a( 2 1i ' ' ' > ^ni^ii • • - ^ n ), so it follows that s T x is also invariant 
under such permutations. Moreover, ( 122]) is divisible by each tjZi + Zj with i < j 
in the unique factorization ring G[zi, ■ ■ ■ ,z n ,t±,--- ,t n ] . The symmetry property 
implies that it is also divisible by tiZj + Z{ with i < j, and since these are coprime 
to Yli<j(tj z i + z j) ^ follows that Z{&\) is divisible by these. Therefore s T x is a 
polynomial in <C[zi, • • • , z n , ti, • • • ,t n ]. 

It remains to be seen that is independent of the In 

A n i <i(^i+^)' 

we regard the numerator and the denominator as both being elements of R\tj\ where 
R = <D[zi, • • • ,z n ,tj(j 7^ i)]. From what we have shown, is a polynomial. We 
claim that both the numerator and denominator have the same degree i — 1 in t{. 
For the denominator, this is clear. For the numerator, the number of — in the top 
row of vertical lattice edge spins is n by the boundary conditions, and it follows 
from Lemma [2] that each successive row has one fewer — . This means that there 
are i — 1 vertices labeled i such that the spin on the edge below it is — , and from 
Table [H it follows that the number of Boltzmann weights equal to Zi(t,i + 1) or tj in 
any particular state is ^ % — 1. The degree of the numerator is thus ^ i — 1 and 
since the degree of the denominator is i — 1, and the quotient is a polynomial, both 
numerator and denominator must have degree % — 1 in t^. Thus the quotient has 
degree zero, and does not involve t^. □ 
We may now conclude the proof of Theorem [5] by showing that = s\. Since 
is independent of ti, we may take all = —1. Now in (1211) the denominator 
becomes fli-cyC 2 ^ — z i)- Since this is skew-symmetric under permutations, the nu- 
merator Z(& v x ) is also skew-symmetric. With tj = —1 any state containing a vertex 

in configuration has Boltzmann weight 0, so we are limited to states omitting 
this configuration. In view of the bijection between states and strict Gelfand-Tsetlin 
patterns, this means that the corresponding Gelfand-Tsetlin pattern T has the prop- 
erty that every entry from any row but the first is equal to one of the two entries 
directly above it. It is easy to see that the weight /i of such a coefficient, described 
by Lemma EJ is a permutation a of the top row of T, that is, of A + p. These weights 
are all distinct since A + p is strongly dominant, i.e. without repeated entries. Since 
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l-h 
3 



Z 3 



Pa{j)+K{ 3 ) 



it is skew-symmetric, its value is sgn(<r) times a constant times \\ z 
To determine the constant, we may take the state whose Gelfand-Tsetlin pattern is 



Ai + pi A 2 + P2 

A 2 + P2 



A n 



A, 



> . 



This has weight Yl z j +Pl an d so 



s \( z li ' ' ' ! z n) 



UiKjiZj-Zi) 



which equals s\(zi, ■ ■ ■ , z n ) by the Weyl character formula. 



6 Tokuyama's theorem 

We recall some definitions from Tokuyama [37] . An entry of a Gelfand-Tsetlin pat- 
tern (not in the top row) is classified as left-leaning if it equals the entry above it 
and to the left. It is right-leaning if it equals the entry above it and to the right. It 
is special if it is neither left- nor right-leaning. Thus in ( 120]) , the 3 in the bottom row 
is left-leaning, the in the second row is right-leaning and the 3 in the middle row 
is special. If X is a Gelfand-Tsetlin pattern, let be the number of left-leaning 
entries. Let <4(T) be the sum of the k-th row of X, and rf„ + i(T) = 0. 

Theorem 6 (Tokuyama) We have 

E f II t^\t + = HO* + tzfrM, ■ ■ ■ ,*n), 

% \fe=l / i<j 

where the sum is over all strict Gelfand-Tsetlin patterns with top row A + p. 

Proof If T corresponds to a state of the Gamma ice with boundary conditions 
determined by A, then we will show that the Boltzmann weight of the state is the 
term on the left-hand side. From Lemma [3] the powers of z are correct. It is easy 
to see that if an entry in the k-th row of X is left leaning (respectively special), and 
that entry is j, then the configuration in the j-column and the k-th row of the ice is 

9. $ 

0«'0 O » '(D 

^ respectively ^ 
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so from Table [U it follows that the powers of t{ are also correct. The statement now 
follows from Theorem [5j □ 



7 The Yang-Baxter equation for Gamma- Gamma 
ice 

We will prove a star-triangle relation that only involves Gamma-Gamma ice. Let 
us think of Gamma ice as being organized into strands of horizontal lattice edges, 
with every Gamma vertex of the strand having the same label i. We may think of 
Gamma-Gamma ice as a tool that switches two strands. The following result states 
that this tool respects the braid relation. We have drawn this picture differently from 
that in Theorem H] since this Yang-Baxter equation involves only horizontal edges, 
while that in Theorem H] involves both horizontal and vertical edges. 



With <7, t, (3, a, p,6 fixed, we may regard these two configurations as ensembles each 
involving three Gamma-Gamma vertices. The Yang-Baxter equation says that they 
have the same partition function. 

Theorem 7 The Yang-Baxter equation is true in the form 



with R = Rrr ■ 

Proof This follows from Theorem [3] since TiYvihj) — const X7Tr(i)7rr(j) _1 , so 




m ky«R(i, k)%R(i, s)z = £ ru, ao&rg, n 



Bp 




7Trr(«, iKirC?, k) = const X7r rr (z, k). 



□ 
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8 More Star- Triangle Relations 



There are further star-triangle relations which go outside the six-vertex model. We 
find that the discussion in Section [I] can be extended the set of Boltzmann weights 
in the eight vertex model that has either aia 2 + b\b 2 — c\c 2 = and d\ = d 2 = 
or a±a 2 + b\b 2 — d\d 2 = and c\ = c 2 = 0. The parameter subgroup will have the 
GL(2, C) x GL(1, C) of Theorem [3] as a subgroup of index two. Let TZ be the set of 
R as in (jSJ) with such weights, where it is assumed a\a 2 + b x b 2 ^ 0. 

Theorem 8 There exists a composition law on TZ* such that if R,T G TZ* , and if 
S = RoT is the composition then [i?, S,TJ =0. This composition law is determined 
by the condition that ir(S) = tt(R)it(T) where it : TZ* — > GL(4, C) is the map defined 
by ljm\ ) if Ci,c 2 are nonzero, and by 



7T(R) = 7T 



ifdi,d 2 are nonzero. 



( ai 



\ d 2 



di \ 



«2 / 



/ 



V d 2 



di \ 



ia 2 
ib 2 



—ibi 
iai 



Here % = y— 1. 

Proof Let us call R £ TZ* of Type C if ci,c 2 are nonzero (so d\ = d 2 = 0) and of 
Type D in the other case. There are four cases to consider. One, where R and T 
are both of type C, is already in Theorem El In the other three cases, we compute 
{R, S,TJ =0 with S as follows. 

If R is of type C and T is of type D then S is of type D with 



a x {S) = a 2 (R)a 1 (T) + b 1 (R)b 1 (T), 

a 2 (S) = ai (R)a 2 (T) + b 2 (R)b 2 (T), 

h(S) = -b 2 (R)a 1 (T) + a 1 (R)b 1 (T), 

b 2 (S) = -6 1 (i?)a 2 (T) + a 2 ( J R)6 2 (T), 

d x (S) = Citi^CT), 

d 2 (S) = c 2 (R)d 1 (T). 
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If R is of type D and T is of type C then S is of type D with 



ai (S) = ai {R)a 2 {T) + b 2 {R)b 2 {T), 

a 2 (S) = a 2 (R)a 1 (T) + b 1 (R)b 1 (T), 

h(S) = b 1 (R)a 1 (T)-a 2 (R)b 2 (T), 

b 2 (S) = b 2 (R)ai(T) — ai(R)bi(T), 

d x (S) = d 1 {R)c 2 {T), 

d 2 (S) = d 2 (R) Cl (T). 

Finally, if R and T are of type D then S is of type C with 

ai(S) = -a2(i2)o 2 (T) + 6i(i2)62(T), 

a 2 (S) = -a 1 (R)a 1 (T)-b 2 (R)b 1 (T), 

h(S) = b 2 (R)a 2 (T) + ai (R)b 2 (T), 

b 2 (S) = b 1 (R)a 1 (T) + a 2 (R)b 1 (T), 

Cl (S) = d l {R)d 2 {T), 

c 2 (S) = d^d^T). 

These computations may be translated into the identity 7r(S) = tt(R) n(T). □ 
We will give some applications of this. The Boltzmann weights for a variety of 
other models are given in Table [2j While Gamma ice is of Type C in the terminology 
of the last proof, we also introduce Delta ice which is of Type D. Delta-Delta ice is of 
Type D and Gamma-Delta and Delta-Gamma ice are of Type C. We will distinguish 
between Gamma ice and Delta ice by using • to represent Gamma ice and o to 
represent Delta ice, and variants of this convention will also distinguish the other 
four types of ice. 

Thus in addition to (jT8l) we have: 



' Zi^i ~\~ Zj 



-Raa(^) tii Zj, tj t 



RrA(Zi, ti, Zj, tj) 



' Zi ~\~ titjZj 



\ Z jtj + Z j 



Zr^ Zj Zj "t/j | Zj 

Ziti ~\~ Z{ Zjtj Z^ti 



Zjtj ~\~ Zi 



Zjti H~ z% 



Zjtj + Zi J 



iti ~\~ Zi ^ 



Zi Zj J 
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Table 2: Boltzmann weights for various types of ice with spectral parameters (zj,tj 
and (zj,tj). (See Table [T]for Gamma and Gamma-Gamma ice.) 



/ Zi- 



R&r( z i, ti, Zj, t J/ 



jtj + Zj \ 



Zyt j I Zj 

\ ZiXi ~t~ ^2 tjtjZi J 

We will denote by r(zj, U) what was previously denoted T(i). We have also A(zj, t 

( Zi 1 \ 



/ 1 \ 

U (U + Y)Zi 
1 Zi 



A(zi,ti 



\ z^U + 1) 1 J 
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Theorem 9 If X,Y e {T, A} then 

iRxYfaUiZj^^faU^Yfatj)} =0. (23) 

Proof In each of the four cases 

-K(R X Y{zi,t u z j ,t j ))'K(Y(z j ,t j )) = Zj(tj + 1) x ir(X(zi,ti)). 

The result then follows from Theorem EJ □ 
Now we turn to generalizations of the Yang-Baxter equation. For every choice 
of z and t and X 6 {T,A}, let V x (z,t) be a two-dimensional vector space with 
basis v x (z,t) and v Y (z,t). Then R XY (zi, £ 1; z 2 , t 2 ) defines an endomorphism of 
V x (z 1 ,t 1 )®V Y (z 2 ,t 2 )by 

R(v a ®v T ) = J2 R *>v®v fl , R = R XY {z 1 ,t 1 ,z 2 ,t 2 ). 
Theorem 10 IfX,Y,Z e {T, A} then we have 

{Rxy( z 1i tl, z 2, h), Rxz( z li tl, z 3, h), Ryz( z 2, h, z Zi = 0. 

Moreover 

{Rxy( z 2i tl, z l,t 2 ), Rxz( z 3i t\i z \i h), Ryz{ z Zi h, z 2i H)\ = 0- 

Proof This follows from Theorem [HJ □ 
We now describe the boundary conditions for Delta ice in the ensemble that 
appears in the second identity in (J3J). The columns are labeled, as with the Gamma 
ice, in decreasing order. However we label the vertices in decreasing row order, so 
the labels of the vertices of the top row are n, and so forth. 

The Delta ice boundary conditions are as follows. We again fix a partition A. On 
the left boundary edges, we put — ; on the right and bottom edges we put +. On 
the top, we put — at every column labeled Aj + n — % (1 ^ i ^ n), that is, for the 
columns labeled with values in A + p. Top edges not labeled by Aj + n — i for any i 
are given spin +. Thus if A = (3, 1, 0), here is the Delta ice. (To indicate that this 
is Delta ice, the vertices are marked o.) 
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Theorem 11 The partition function is 

Z(ef)(z!, ■■■ ,z n ;t h --- , t n ) = Y\(tjZj + z i )s x (z 1 , 



Z n )- 



i<j 



Proof This is proved analogously to Theorem [5J using the case X 
Theorem [91 We leave the details of the proof to the reader. 
Theorem [9] may be used to show that 

Y[(t jZj + Zi )z{& v x ) = z(e$) Ufa + zi) 



Y 



A of 

□ 



(24) 



i<j 



directly without invoking Theorems [5] and [TTJ This fact is closely related to State- 
ment B in Brubaker, Bump and Friedberg [3], and the following argument may be 
used to give an alternative proof of that result in the special case where the degree 
(denoted n in [3]) equals 1. 

Begin with an element x of & T Xl say (for example with A = (3, 1, 0)): 



1 
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(The unlabeled edges can be filled in arbitrarily.) We wish to transform this into an 
element of an ensemble that has a row of Delta ice so that we may use the mixed 
star-triangle relation. We simply change the signs of all the entries on the edges in 
the 3 row: 

5 4 3 2 1 

Q $ $ Q $ Q 

©fofofofofofe 


(±) (±) o @ @ @ 

£D> £D> £D> £D> £D> £Di 
\D \D \D \D \D \D 

Let x' be this element of the mixed ensemble & . We observe that the Boltzmann 
weights satisfy w(x) = w(x'). Indeed, in the bottom row only the following types of 
Gamma ice can appear: 



These change to: 



Gamma 
Ice 


0«'0 


(D? 2 






1 


l 






Delta 
Ice 


e^e 




$ 




i 


i + 





Observe that the weights are unchanged. Note that this would not work in any 
row but the last because it is essential that there be no — on the bottom edge spins. 
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Now we add a Gamma-Delta R-vertex. 

8 8 8 8 8 8 

©fofofofofofe 


^>?ofofofofofe 

(T\ (T\ (T\ (T\ (T\ 
\D \D \D vf vf vf 

If 6" is this ensemble, we claim that = (t 3 z 3 + z 2 )^(6') = (hz 3 + z 2 )Z(&\). 

Indeed, from Table [2l the values of a and b must be +, — respectively and so the 
value of the R-vertex is t 3 z 3 + z 2 for every element of the ensemble. Now using the 
star-triangle relation, we obtain Z{&") = Z(&"') where &"' is the ensemble: 

5 4 3 2 1 




8 8 8 8 & 8 



Here we must have c, d = +, -and so {t 3 z 3 + z 2 ) Z (& T X ) = Z{&") = (t 2 z 3 + z 2 )Z(e {iv) ) 
where &^ is the ensemble: 



1 







8 8 8 8 8 8 

a a a a a a 

©^o^ofofofofe 


8 8 8 8 8 8 



29 



We repeat the process, first moving the Delta layer up to the top, then introducing 
another Delta layer at the bottom, etc., until we have the ensemble <5^, obtain- 
ing (EU). 



9 Yang-Baxter Systems 



The results of this section are further applications of Theorem [HJ 

An important property of the R-matrices Rxy{zi,ti, Zj,tj) is that they are pro- 
jectively triangular. That is, 



RxY( z i,U, Zj,tj) 1 — CxY(Zi,ti, Zj,tj)P Ryx( z j,tj, Zi,ti) P 

where cxy( z i,U, z j^j) is a scalar and 

/ 1 



(25) 



V 



1 7 



The constant cxy m &y be eliminated by multiplying Rxy by a suitable scalar - for 
example in the case X = Y = T if R[* v (zi,ti, Zj,tj) = (zjU + Zij^Rrri^ijU, z j,tj) 
then i?p r satisfies ( |25l) without the cxy, at the cost of introducing denominators. 

Yang-Baxter systems occur with varying degrees of generality in connection with 
different problems. One type occurs in the work of Vladimirov [39] on quantum 
doubles; another type occurs in Hlavaty jH] on quantized braided groups. The most 
general formulation [T5] , [13] involves four types of matrices which correspond to our 
Rxy, X,Y e {r, A}. 

The axioms for a parametrized (or "colored") Yang-Baxter system in the most 
general definition require four types of matrices, A, B, C, D, depending on parameters 
z\ and z 2 and subject to the properties 



[A A, Aj 
{A, C, Cj 

A, C, B$ 



0, 
0, 
0, 
0, 



{D,D,D} 
{D,B,B\ 
D,C\C$ 
D,B,C$ 



0. 
0. 
0. 
0. 



(26) 



where we now denote 

{X, Y, Zj = X 12 (zi, z 2 )Y 13 (z h z 3 )Z 2 3(z 2 , z 3 ) - Z 23 (z 2 , z 3 )Y 13 (zx, z 3 )X 12 (z h z 2 ) 
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and X*(zx,z 2 ) = PX(z 2 , Z\)P. We have two spectral parameters z and t, so we 
interpret 

X i (z 1 ,t 1 ,z 2 ,t 2 ) = PX(z 2 ,t 2 ,z 1 ,t 1 )P. 

Theorem 12 Let X,Y e {r,A}. Then 

A = R XX , C = B$ = R XY , D = R YY X 
is a Yang-Baxter system satisfying (E6jj. 

Proof We leave the verification to the reader. □ 
Note that by projective triangularity we may replace B by Ryx~, which is a 
scalar multiple of Rxy X ■ Thus if X = T, Y = A we have the Yang-Baxter system 

A = Rrr, B = R&r 1 , C = i?rA, D = -Raa^i 

which uses each of the four braided ice types in Table [2] exactly once. It is probably 
most interesting to take 1^7, but worth noting that we can also make a Yang- 
Baxter system with i?rr (or Raa) playing all four roles. And we also obtain a Yang- 
Baxter system as follows by interchanging the Z{ (but not the £j) in the spectral 
parameters. 

Theorem 13 Another set of four Yang-Baxter systems may be obtained by taking 
A = Rxx, C = B* = Rxy, D = Ry Y , 

where 

RxY{Zi,ti, z 2 , t 2 ) = Rxy(z 2 , ti, Zi,t 2 ). 
Proof We leave this to the reader. □ 
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